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Abstract. The well-known theorems of Stieltjes, Hamburger and Hausdorff establish con- 
ditions on infinite sequences of real numbers to be moment sequences. Further, works by 
Caratheodory, Schur and Nevanlinna connect moment problems to problems in function 
^ , theory and functions belonging to various spaces. In many problems associated with re- 

■^^ 1 alization of a signal or an image, data may be corrupt or missing. Reconstruction of a 

function from moment sequences with missing terms is an interesting problem leading to 

advances in image and/or signal reconstruction. It is easy to show that a subsequence of a 

moment sequence may not be a moment sequence. Conditions are obtained to show how 

rigid the space of sub-moment sequences is and necessary and sufficient conditions for a 

^^ ' sequence to be a sub-moment sequence are estabhshed. A deep connection between the 

^ ' sub-moment measures and the moment measures is derived and the determinacy of the mo- 

V | _ | J , ment and sub-moment problems are related. This problem is further related to completion 

_5 . of positive Hankel matrices. 

-)— > 

B- 

1. Introduction 

CN , While numerous applications of moment problems with a complete set of moments have 

1^ ' been been identified, they are mostly limited to theoretical observations. For practical im- 

rjL , plementation of moment problems, it is vital to be able to deal with missing moment data 

QQ ' since data obtained from physical experiments and phenomena are often corrupt or incom- 

l/^ , plete. We call a positive subsequence of a moment sequence a sub-moment sequence and 

-Y-\ ' the moment problem of sub-moment sequence a sub-moment problem. In this paper we 

(^ , give a complete characterization of the sub-moment problem and develop a deep connec- 

CsJ ' tion between the measures arising from the sub-moment problem and the original moment 

problem. Further, we relate this to completion of positive Hankel matrices and arrive at a 
necessary and sufficient condition for the completion of such matrices. Recent advances 
vj , in image reconstruction from sparse MRI data |6|, partial data transmission and sparse 

rS ' signals, make this problem timely and of broad interest. 

C^ , In his memoir "Recherches sur les Fractions Continues" [SO] from 1894-95, Stieltjes 

introduced the "Problem of Moments" on [0,oo): find a bounded non-decreasing function 
o-(u) in the interval [0,oo) such that its "moments", s„, are given as 

r*co 

I M'Vcr(M), «elNo = {0,l,2,---) 
Jo 

Stieltjes took the terminology "Problem of Moments" from Mechanics. He often used 

mechanical concepts of mass, density, stability, etc., in solving analytical problems. We 

r 

can consider dcr(u) to be mass distributed over the interval [u,u + du] so that I dcr(t) 

Jo 

gives the mass distributed over [0,m]. Accordingly, I udcr{u) and I m £/cr(M) respec- 

Jo Jo 

lively represent the first moment and the second moment, also called moment of inertia. 
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with respect to of the total mass I dcr(u) distributed over the interval [0, oo). Stieltjes 

Jo 

/-•oo 

generalizes this concept to call I u"da-(u) the n-th moment of the given mass distribution 

Jo 
characterized by the function o-(u). Stieltjes proves several results concerning solubility of 

the problem of moments and uniqueness of the solution. He shows that for a sequence {s„} 

the moment problem is solvable if and only if the determinants of the matrices 
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I S„+i ■■■ S2n \ L^H+l Sn+2 '" ■Sln+l. 

are non-negative. There have been several sufficiency conditions established for unique- 
ness of the solution to the moment problem. While Stieltjes had mechanics in mind, Cheby- 
shev and Markov had explicitly probability theory as aims of moment problems. 

Hamburger continued the work of Stieltjes in his series of papers "Uber eine Erweiterung 
des Stieltjes Momentenproblems lfT6l ". from 1920-21. He treated the moment problems 
as a theory of its own. He extended the Stieltjes' moment problem to the whole real axis. 
The Hamburger Moment Problem is stated as: find a non-decreasing function o-{u) in the 
interval (-00,00) such that its moments, s„, are given by 



Jk 



u"dcr{u), n e INq. 



He established that a necessary and suflicient condition for the moment solution to exist 
is the positivity of moment sequence. Later, mathematicians including R. Nevanlinna, M. 
Riesz, T. Carleman, F. Hausdorff, M. Stone, and C. Caratheodory further studied several 
moment problems. 

A moment problem that is not well-posed may not be solvable. Criteria for existence of 
solution is studied to a great extent for every type of moment problem. Once the solubility 
criteria are determined, then comes the question of uniqueness of the solution since it is 
often the case that there is more than one solution. Two solutions of a moment problem are 
not considered to be distinct if their difference is a constant at all the points of continuity 
of the difference. Such a moment problem with no more than one distinct solution is called 
determinate, and indeterminate otherwise. 

Hausdorff 's moment problem deals with finding a function a-{u) with support on the 
closed unit interval [0, 1] such that a given sequence {s„] satisfies 

-1 



Sn - I u"o-(u), n e INq. 
Jo 



Hausdorff establishes that there is a unique function cr(u) in [0, 1] if and only if the se- 
quence {sk} is completely monotonic, that is, it satisfies 

(-l)"A"sk>0, foralln,fe>0, 

where A is the difference operator defined as 



(=0 ^ ' ' 



This condition for existence and uniqueness holds for the moment problem in any finite 
interval. In fact, the Hausdorff moment problem is often stated on the interval [-1, 1] just 
for the sake of simplicity. 
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The problem of moments stands as a very important problem in analysis up to the 
present day. Results on moment problems have numerous applications in the areas of 
extremal problems, optimization theory and limit theorems in probability theory. Rather 
than the usual infinite moment sequence problem, the truncated moment problem is more 
appUcable in mathematical and physical sciences. Further, multivariable moment problem 
is more applicable than a single variable one. Multivariate extension of moment problem 
is also largely studied (see [11] and |26| for instance). 

Moments of a geometric object conveys a lot of its topological information such as how 
round it is, where it is located, what direction it is tapered, where its mass is centered, 
etc. For simple objects like an ellipsoid, just a few terms in the moment sequence are 
enough to identify it. However, only infinitely many terms can uniquely identify more 
complicated shapes like a polygon or a quadrature domain. The authors in fZSl give an 
example of identifying a polygon from its given moments. They prove that the vertices 
of a nondegenerate, simply connected n-gon in the plane are uniquely determined by its 
moments up to order 2«- 3. Theory of moments can also be applied to some modeling 
problems. An example of speech modeling can be found in [9J. Additional examples of 
data reconstruction with moments can be found in lITOl . llT4l . ifTsl . llT9l and ll29l . 

Various generalizations and extensions of the problems of moments have been studied 
to deal with several applications. Some generalizations include replacing the sequence of 
functions {u"} with a more general sequence /„(m) such as {e'""}, and replacing integrals 
with more general functionals in abstract spaces. Lasser [23 1, for example, studies moment 
problem by replacing {u"} with polynomials. However, we will not consider these gener- 
alizations in this paper The important variation of the moment problem we look at is the 
moment problem with some of the moments missing. Moment data obtained from various 
experiments are often corrupt and incomplete. Such reconstruction of moments can be then 
applied in problems associated with realization of a signal or an image. Reconstruction of 
a function from moment sequences with missing terms is an interesting problem leading to 
several advances in image and/or signal reconstruction. In this paper we will work towards 
finding techniques for reconstructing the missing moment data. 

Several authors have studied this problem by replacing the missing moments with an 
appropriate value, leading to a perturbation of the moment sequences followed by the 
question of sensitivity of the corresponding orthogonal polynomials. Because of necessary 
requirements for a sequence to be moments, this approach turns out to be very restrictive in 
choosing and replacing the missing moments. Thus, we look for ideas to extract a positive 
subsequence of a positive sequence and look for its moment solution. 

We primarily have two questions to answer: 

Ql. How do we characterize positive subsequences of a positive sequence? 

Q2. How do we characterize the corresponding non-decreasing measures arising from 

the sub-moment data and what, if any, is the relationship of these measures to the 

original measures? 

In this paper we establish several methods for extracting positive subsequences and look 
at solutions to the sub-moment problems. Unless stated otherwise, by moment problem we 
shall refer to the Hamburger moment problem throughout the paper In the first section 
of this paper, we will thoroughly discuss some classical results and properties of positive 
sequences and corresponding polynomials. Then we will survey some of the prominent 
works in the field of moment problem in the next couple of sections. The results described 
here, mostly without proof and with reference to literature, are those that will play an es- 
sential role in the sub-moment problem. The second part of this paper is dedicated to the 
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Study of the sub-moment sequences, sub-moment problems and positive Hankel matrix 
completion problem. We shall establish necessary and sufficient conditions for a subse- 
quence of a moment sequence to be positive. We will also find conditions for determinacy 
of sub-moment problems. These results give a complete answer to the problem of Hankel 
matrix completion. 

2. Positive Sequences 

Positive sequences play an important role in the theory of orthogonal polynomials, pos- 
itive functionals and moment problems. This section is concerned with an overview of 
some characteristics of positive sequences. In this section, we will describe many of the 
classical theorems and results that we will need in the ensuing sections of this paper In 
most cases the proofs are not given and the reader is referred to one of the several paper or 
classical monographs on the theory. In a few instances, for later reference and for readers' 
convenience, brief proofs are provided. 

Definition 2.1. An infinite sequence {s^] is positive if the quadratic form 



/ 1 ^i^j^ 
iJ=Q 



'+.I 



is positive for any m and xq,x\,--- ,x„, e R. This is equivalent to saying that the determi- 
nant 



(1) 



£>,„ = 



^0 






^2w 



is non-negative for m e Mq. 

Proposition 2.2. Term-wise sum and product of two positive sequences is positive. 

Proof. The claim that the term-wise sum is positive follows readily from the definition. 
The claim that the term-wise product is positive follows from the Schur product theorem 
IfTSl Theorem 7.5.31. n 

Given a sequence {s^}, define a linear functional A such that 

A(P(m)) = P0S0-\-piSi + ■ ■ --^PnSn 

where 

P{u) — pQ + piu H v p„u^ . 

Definition 2.3. A linear functional A is called positive if for any polynomial f : IR — > R, 
P{u) > and P(u) ^ 0, it follows that A(P(u)) > 0. 

Tlieorem 2.4. For any sequence {s^}, the functional A defined above is positive if and only 
if{sk} is positive. 

Proof. Suppose {sk} is positive. Any polynomial P(u) - po + piu H — -^ Pnu" > and 
P{u) t- can be written in the form 



P(m) = [A(m)]2h-[B(m)]2 
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where A(u),B(u) are polynomials with real coefficients. Then we have 

m m 

A(P(u)) = 2_^ xiXjSi+j + 2_] ytyjSi+j > 0. 

i,j=0 i,j=0 

Now suppose A is positive. Then 

m 

J]xiXjSi+j = A([P(u)f), 

ij=0 

where P(u) = xq + xim h h x„,m"'. Hence {sk} is positive. □ 

Given a positive sequence {.s^), we can construct a set of polynomials {P„{u)} such that 
P„(u) is of degree n with positive leading coefficient and they are orthonormal with respect 



to the functional A. These polynomials are explicit 



(2) 



where 



Pn(u) ■■ 



1 



VA,-i£>« 



Si 
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y given by the formula 



A, = 



so 
si 
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Sn-l 
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Sn 
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S2n S2n-l 
,,n—l ,,« 



andD_ 



It is well known that this orthonormal set of polynomials {P„(m)) satisfy the three term 
recurrence relation 



(3) 
where 

and 



uPn{u)-bnPn+l{u) + anPn{u) + b„-iPn-l{u), n- 1,2,3, 

yiDk-iDk+i 



a„ - A(u[P„(u)]n , and b„ - 
Po(m) = 1, andPi(M) = 



u — Qq 



Together with the orthonormal polynomials {P„(m)), the polynomials of the second kind 
{Qn{u)] are also generated from the same three-term recurrence relation (|3]l with initial 

1 
conditions Qq{u) - Q,Qi{u) - — . {Pn{u)) and {Qn{u)} are linearly independent solutions 

_ bo 

to ^. In fact, Qniu) is a polynomial of degree « - 1 and can be expressed explicitly as 



(4) 



Qniu) = A 



Pn{u)-Pn{v) 



Orthonormal sets of polynomials {P„(m)) and {Qn{u)] play an important role in the the- 
ory of moment problems. See [13] for further details on the polynomials and their recur- 
rence relation. Several conditions on these polynomials have been established to determine 
the existence and uniqueness of solutions of moment problems. We will elaborate more on 
the concepts and usage of these polynomials in the later sections. 
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3. Classical Moment Problems & Determinacy 

In this section we review some of the key resuhs in the Hamburger's moment prob- 
lem, also called power moment problem, which is stated as: given an infinite sequence of 
numbers {sk}, find a non-decreasing function o-(u) such that 



(5) Sk = I u cr(u) 



= I u o-(u) 
JR 



This problem was formulated and solved by Hamburger Later works by mathematicians 
including Fischer, Akheizer and Krein establish conditions on the spectrum of the moment 
solutions. Some of these results are summarized in the theorems below. The reader is 
referred ||T] and ||281 for an extensive account of the principle results. 

Theorem 3.1 (Hamburger). There is a non-decreasing function cr(u) with an infinite num- 
ber of points of increase such that equation Q holds if and only if{sk} is positive. 

Theorem 3.2. For a non-decreasing function cr(u) with precisely A:-H 1 distinct points of 
increase such that equation (|5]l holds, it is necessary and sufficient that 

Do,Du--,Dk>0, DM^Dk+2 = ■■■ = 0, 

where Di is the determinant of the Hankel matrix defined in ([T]). 

Foraproof of the Hamburger Theorem, refer to fll]. The proof uses the Helly's Theorem 
extensively. In some instances, such as in fST] and pTl, the properties of the determinants 
Dm have been studied to describe the moment solution. 

The orthonormal set of polynomials {Pn{u)] and {Qn{u)] are used to establish the con- 
ditions on determinacy of the moment problems. So before moving on to results on deter- 
minacy we will define the circle Kn{X) as follows. 

Theorem 3.3 (Helhnger). Let AeG be fixed with 3/1 > (ViA < Q) and let t e R. Then 

,.. ,, , Qn(^)-TQ„-i(A) Q„iA,T) 

(o) w„{A,t) — — 

Pn(A)-TP„-i(A) P„(A,t) 

describes a circular contour K„(A) in the open half plane 3w > (3w < 0). The center of 
this circle is at the point 



^7^ Qn{A)P„.i{A)-Q„.,{A)P„{A) 

P„{A)P„-i(A)-Pn-i{A)PjX) 
with radius 

m - ' 



n-l 

\A-J\Y^\Pk(At 
Furthermore, the equation of the circle can be written in the form 

_— n-l 

(9) '^^-Y,\''PkiX) + QkiAt^Q 

Proof Refer to ETl. □ 

Note that /r„+i(/l) c K^iA). We call the limiting circle KoaiA)- It can be shown that the 
quantity 

) 



(10) w{A) 

In u-A 



f dcr(u) 
Jth u-A 
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falls on the circumference of Kco(A), |lj p. 30-34]. 

Now we look at some results concerning determinacy of the Hamburger moment prob- 
lem. Most of the determinacy theorems we talk about here deal with the nature of the 
limiting circle Kca(A). It is useful to present the following theorem and its proof on deter- 
minacy since we will be using it for the case of sub-moment problem in a later section. 

Theorem 3.4. If {sk} is a positive sequence such that its corresponding limiting circles 
Kco{A) is a point, then the Hamburger moment problem of sequence {s^} is determinate. 

Proof. Suppose the moment problem has two solutions cri{u) and a-i{u). Consider the 
functions 

,,,s ., . r do-i(u) ,„,. r da-2(u) 

(11) /i(w)= —, and fiiu)^ -, 

J^ u-A J^ u-A 

for some Ae<C such that V>A + 0. Now for a fixed A, each of /i(/l) and /2(/l) are in Kca{X) 
by Theorem 2.2.4 in |1 1. But \iKaa{A^ is a point, then/i(/i) and/2(/l) coincide. Putting 

a(u) — cri(M) - cr2(u), 

we have 



Jr 



da(u) ^ 



JR u-A 
By the Stieltjes and Perron inversion formula, 

a(u -0) + a(u + 0) 

= c for some constant c. 

2 

Therefore criiu) and crjiu) represent the same solution of the moment problem. n 

The converse of this theorem is also true. That is, if the Umiting circle K^{A) is not a 
point, then the Hamburger moment problem is indeterminate. Carlemanf? | proved several 
striking results regarding determinacy of the Hamburger moment problem. We mention a 
couple of them here. 

Theorem 3.5 (Carleman). If for a positive sequence {si^}, 

= oo, 

where bn are constants as given in the three term recurrence relation ([3]l, then the moment 
problem is determinate. 

Theorem 3.6 (Carleman). If for a positive sequence {sk], 

oo ^ J 

n=l 

then the moment problem is determinate. 

4. Function Theory and the Moment Problems 

For a measure criu) which arises as a moment solution to a positive sequence, we can 
define the V space in the usual sense: 

(12) L^ = < / : R — > C : / /i cr - measurable and I \f{u)\''dcr(u') <oo\ . 
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The space of functions L^ with the norm 

(13) \\f(u)\\^,^=lj \f(urdCT(u)Y 

is complete and is a Hilbert Space when p = 2 where inner product defined as 

(14) (f(u),g(u))^^ f f(u)'^dcr(u) 

Jn 

for any functions f,g € L^. 

Use of the space L^ to study the properties such as determinacy of the solution(s) <t(u) 
for the moment problem is very prevalent. This is done mostly in connection to the poly- 
nomials generated by the corresponding sequence {sk}. Several theorems by M. Riesz 
characterize solutions of a moment problem and its determinacy ll27l . Below is a very 
useful consequence of these theorems. 

Corollary 4.1. Ifcr is the solution of a determinate moment problem, then the set of all the 
polynomials is dense in L^. 

The Nevanlinna formula relates all the solutions of an indeterminate moment problem to 
the functions in A^-class. This formula establishes a one-to-one correspondence between 
the aggregate of all the solutions cr of an indeterminate moment problem and the aggregate 
of all the functions ^ in class N, which we will define below, augmented by the constant 
at oo. Applying the Stieltjes-Perron inversion formula to the NevanUnna formula one can 
determine cr in terms of ^. 

Definition 4.2. Nevanlinna matrix is any matrix of the form 

a{z) b(z) 
c(z) d(z) 

where the elements a(z),b(z),c{z) are entire transcendental functions, and the following 
two conditions are satisfied: 

(1) a(z)d(z)-biz)c(z)=l 

(2) For any fixed real number t e IR the function 

Ta(z)-c(z) 
Tb(z) - d(z) 

is regular in each open half plane and satisfies 

3w(z) 



9z 



>0. 
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Now, let us consider the four families of polynomials given as 

n-l 
Aniz)^zY,QkiQ)Qk{z) 



k=0 



Bn(z)^-l+zY,Qk(0)Pk(z) 
k=Q 
n-l 

Cn(z)^l+zY,PkiO)Qk(z) 



k=0 
n-l 

Dn(z)^zJ]Pk(0)Pk(z). 
k=0 

If the moment problem is indeterminate, then the sums 

CX) CO 

Y,\Pk(z)\\ Y,\Qk(z)\^ 

k=Q k=Q 

are bounded and converge uniformly to respective limits on compact subsets of C. By 
the Cauchy-Bunyakovskii inequality, the polynomials A„(z),B„(z),C„{z),Dn{z) converge 
uniformly to the entire transcendental functions 

cx> 

Aiz)=zY,Qk(0)Qk(z) 

k=0 

CO 

B(z) = -l+zY,Qk(0)Pk(z) 

k=0 

CO 

Ciz)^l+zJ]Pk(0)Qk(z) 

k=Q 

CO 

D(z) ^zY,Pk(0)Pk(z) 

k=Q 

on compact subsets of C. 

It has already been shown that the limit functions A(z),B(z),C(z),D(z) are entire tran- 
scendental. It can also be shown that the functions A(z), B{z), C(z),D(z) satisfy the require- 
ments for elements of a Nevanlinna matrix. Thus the matrix 

'A(z) B{z) 
D(z) D(z) 

is a Nevanlinna matrix. 

Before we state the Nevanlinna formula, we define the Nevanlinna class, N-class, of 
functions. 

Definition 4.3 (N-class). The Nevanlinna class (Class N) consists of all holomorphic func- 
tions w — f(z) on the upper half plane such that Viw > Q. By the Riesz-Herglotz integral 
representation formula, a function f & N may be expressed as 

f{z)^liz + v+ I —dT{u) 

J u-z 

where //, v e R,// > and t(u) is a non-decreasing function with bounded variation. 
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Theorem 4.4. If3z>0, for an indeterminate moment problem, 

■do-(u) Aiz)(piz)-Ciz) 



I 



u-z 



B{z)<p(z)-D{zy 



where A(z),B(z),C(z),D(z) form a Nevanlinna matrix, cr is a solution to the concerned 
moment problem, and 4> € N. 

This theorem leads to the question of relation between the moment solutions of subse- 
quences to the functions in the N-class. It is worthwhile to look for the moment solutions 
of subsequences in the same N-class. 



5. Incomplete and Sparse Moment Sequence 

There are two main methods for dealing with missing moments: (i) Perturbation/modfication, 
and (ii) Subsequences. Perturbation or modification method is done by replacing the miss- 
ing data with an appropriate value, whereas the method of subsequences is based on simply 
removing the missing moments. In some applications, people have tried replacing the lost 
moments with zero as well. Because of Hamburger's theorem, in both cases the modified 
data has to be positive to guarantee a moment solution. Let us look at an example of a 
positive sequence and some of its perturbations and subsequences. 

The sequence given by Sk = - — - is a positive sequence since, for every k, determinant 
of Hankel matrix is simply given as 



Dk 



1 
2 
1 
3 
1 
4 



1 



k+1 
1 

1+2 



1 



k+1 k+2 



2k+l 



n 

1=0 



(2/ +1)1 



2/ 



which is positive. So the sequence {s^} is positive and has a moment solution. 
By setting s^ = 0, the second Hankel matrix has determinant 



D2^ 



1 i 

2 



<0. 



This shows how easily we lose positivity of the sequence if perturbations of modifica- 
tions is not done properly. Thus, perturbation is a very restrictive method because of the 
requirement of positivity. Due to this restriction and involvement of extensive matrix com- 
putations, the study of sensitivity of corresponding polynomials and moment solutions 
becomes computationally expensive, if not intractable. In this chapter we will mostly deal 
with the method of subsequences. 
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The subsequence given by s^j - — e {sk} is not a positive sequence since the determi- 



nant of its second Hankel matrix is 



Di 



1 1 



<0 



2 6 
and hence has no solution to its moment problem. This example shows that an arbitrary 
subsequence of a positive sequence is not necessarily positive. Thus we need to look for 
ideas to appropriately remove moments while preserving positivity of subsequences. 

5.1. Perturbation of a Positive Sequence. In this section we will make a few notes re- 
garding perturbations of terms in a positive sequence. Modification of a term in a mo- 
ment sequence very significantly aff'ects its moment solution and orthogonal polynomial it 
generates. Several people have studied the sensitivity of orthogonal polynomials to mod- 
ifications of the moments. Gautschi [12J . for example, studies sensitivity of orthogonal 
polynomials to perturbations in moments. Reader is advised to refer to U and ll22l for 
some more examples on modified moments. 

The fact that ^o > kii> ' = 1,2,3,- ■ ■ , suggests that we can have a wide range of choices 
for So keeping the rest of terms constant and maintaining positivity. However, we are not 
much interested in modifying sq since we can always scale this term. The theorem below 
shows how rigid positive sequences are with respect to a single term 5„, n ^ 0. 

Theorem 5.1. Let {sk}k>Q be a positive sequence. If 

{Sk if ki^n 
T if fc = n 

is positive and s^ < s^, then t — s„. 

Proof. Consider an arbitrary system of number xo,jci,X2,--- ,x„,. Then we can write the 
sum 

m mm 

(15) 2j ^i^j^i+j — 2j ■"■'■"■j^n + 2j ^i^j^i+i — ^• 
,j=0 ,j=0 i,j=0 

i+j=n '+j^« 

Similarly, the positivity of {Sa:)a:>o gives 

m m m 

(16) 2_j ^rXjSi+j = 2_j XiXjT+ 2_] XiXjSi+j > 0. 

i+j=n i+ji^n 

Combining (fTsT i and ( fTSI l, and using Sk < s^, 

m m m 

/ ^ XiXjT— y ^ XjXjSfj + y ^ XiXjSi+j > 0, 





,j=0 ij=0 ij=0 
i+j=n i+j=n 


which implies 




(17) 


m m 

^ XiXjSi+j > [s„ - t] ^ 

i,;=0 i,j=0 
i+i=n 
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Since {sk}k>o is a positive sequence the left side of ( fTTb should be positive for any 
system of numbers xq,x\,X2,--- ,x,„- Now, let s„ > r and choose xo,xi,X2,--- ,Xm-l- Then 

m m 

y XjXj = (« + 1). Then ( fTTl i gives y xiXjSi+j > (n+ l)(s„ - r) > 0. But this is not true 

i,j=0 i.j=0 

i+j=n 

since, scaUng the sequence {sk} to have the the first term as ^o - — — 7^ — and using 



the fact that st < sq for all /, we have 

m m 

2_] XiXjSi+j < 2_] XiXjSo = {m + 1) 

/,;=0 i,j=0 

Therefore s„ < t. 

Similarly using (fTSl l and (fTSI l again we can obtain 



2(m+l)2 



^0 = ;; <{n + \)(sn-T). 



(18) 



/ J X-iXjSi+j > [t— Sf,\ y ^ XiXj. 



i,j=0 



ij=0 
l+j=n 



Then using the same argument as above we see that s„ > r. This completes the proof, n 

It is noteworthy that if we perturb more than one term then we have more flexibility. The 
moment sequence tends to get less rigid as we replace more terms. If the moments si and 
S2 are replaced by xi and X2 respectively, does imposing positivity imply si - xi,S2 - X2^ 
Suppose Sp.,i= 1 , 2, . . . , n are replaced by Tp. . Then we can write the sum 



(19) 



/ J Xi^j^i+j ~ Z_i Z_i ^'^J^Pq + /_i ^iXj^i+j > 0. 



i,j=0 q=l i,j=0 

•+J=Pi, 

Similarly, the positivity of {sk]k>o gives 



i,j=0 
i+ji=Pl,...,p„ 



(20) 



III It III III 



i,j=0 



=1 ij=0 

i+i=Pq 



i.j=0 

l + i*P\,-;Pl, 



Combining (fT9T l and (|20] |. and using s^ < s^ gives 



(21) 



2 

q=l 



^Pq Z_i ^'^} 
ij=0 

'+j=Pq 



7 i ^i^j^i+j 



i,j=0 



where ^p^^ - Sp^ - Tp^. Now setting Xi - 1 for i -Q,l,...,m from (ISTT i we have 

n 

(22) Yjyp,(p, + l)]<(m+lfso. 

q=\ 

Although very restrictive and not much revealing, perturbation can be useful in some 
cases. An example of an application of modified moments to harmonic solids is given 
in lis). An important application of modification in sq is that it can lead to a determinate 
solution. It was proved in Stieltjes' 1894-Memoir |30| that a determinate moment solution 
can be obtained from an indeterminate one with a modification in sq 
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The moment sequence s„ - q 2 gives the Stiehjes-Wigert polynomials P„(u;q), 
which are orthogonal in a log-normal distribution, known to be indeterminate. The modi- 
fied sequence {S/t)/t>o defined as 

1 

and 5„ - s„ for all n > 1 has a determinate moment solution it(m) given by 

a-{u) ^2_jCu6u, 

u€U 

where U is the zero set of the reproducing kernel 



K(0,w) = J]p„(0)P„(w) 



n=0 

and 

1 

Cu T ,ue U. 

5.2. Positive Subsequences. As was seen in an above example, an arbitrary subsequence 
of a positive sequence is not necessarily positive. We will develop some methods to extract 
positive subsequences from a positive sequence. 

First we give a simple construction of a positive sequence by the following theorem and 
proposition, which can then be used to generate positive subsequences. 

Theorem 5.2. Let fk '■ M. ^> R be a sequence of functions such that 

(23) f{u)fj{u)^f^j{u). 

Then the sequence {sk] given as 

Sk = I fk(u)da-(u) 

is positive for any non-decreasing function cr(u) with an infinite number of points of in- 
crease. 

Proof. For any finite set of numbers xq, xi , ■ ■ ■ , x„, we have 



tn in n 

^ XiXjSi+j = ^ x,x^' I fi+j(u)do-(u 

i,j=0 i,j=Q "^ 

m „ 

= ^ XiXj I f(u)fj(u)do-(u 



) 
i,j=0 

I Y^XiXjfi(u)fj(u)do-(u) 
^ i.j=0 



/m 
Y^Xifi(u) 

.1=0 



do-{u) > 0. 
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m 

Note that the equahty sign holds if the function ^ Xifi(u) is zero at every point where 

/=o 

m 

o-(u) increases. Since cr(u) increases at infinitely many points then the function y Xifi(u) 

/=o 
is identically zero. Therefore we have the equality sign only if jc, = for all /. n 

Proposition 5.3. The sequence fk(u) satisfies the condition fi(u)fj(u) — fi+j(u)for i,j e Mq 
if and only if 

fk(u)^(fl(u)f. 

Proof. Clearly the conclusion holds for k = I. Assume f„{u) - ifi{u))" holds. Then 
fn+\{u) - fn(u)fi(u) - (/i (m))"/i (m) = (/i(m))"''"', inductively. Proof of the converse is triv- 
ial. D 



Note that Theorem 15.21 and Proposition 15.31 allow many constructions of positive se- 
quences. For example, letting fk(u) - u we obtain precisely the classical power moment 
sequence for any non-decreasing measure a-(u). Similarly, letting fk(u) to be a for some 
constant a or ((j)(u)) for any cr-measurable function gives us more positive sequences 
corresponding to a non-decreasing measure cr{u). 

Now we are interested in extracting positive subsequences from a moment sequence. 
Using Theorem 15.21 and Proposition 15.31 we can construct several positive subsequences 
from a moment sequence. For example, let [sk] be a moment sequence constructed from a 
non-decreasing measure o-{u) defined as 

Sk^ j ((P(u)fdo-(u) 

for some function (p{u). Then for a fixed { e INo the sequence 

Skf-.^ I ((0(m))^) do-(u) = I {<p{u)f'^do-{u) = Ski 

is a positive subsequence of {.s^^). But how can we characterize all the positive subsequences 
we can extract from a moment sequence? Given a positive sequence [sk], the problem of 
identifying all the positive subsequences requires finding all the sequences {(k) £ INo such 
that the sequence given as 

Sk - Sk+tt 

is positive. 

Theorem 5.4. Let {sk} be a positive sequence. Then the sequence Sk = ■Sjt+4 is a positive 
sequence if{u *) is a positive sequence for all m € R. 

Proof. By Hamburger's theorem {sk} is a positive sequence is equivalent to the fact that 
there is a non-decreasing function a{u) such that 

Sk - I u^do-{u), k € INo- 
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Forxo,xi,--- ,Xm e K, we have 



i,j=0 i,j=0 

m r. 

i,j=0 -'k 



i,j=0 



(24) 



Vuj=o 



dcr(u) 
dcr(u). 



But the integrand in JTM can be written as 



(25) y u '+J(xiu')(xjU-') = \xo x\u 



««'"] 






m''"' m''""+' 






/2ni 



XO 
XlU 



Thus, the integral (|24] | is positive if the Hankel matrix for the sequence {m *) is positive for 
all M e R. This is equivalent to the sequence {m *) being positive for all m e R. n 

However, the following example due to David Kimsey (private correspondence) shows 
that the converse of this theorem does not hold in general. We gratefully acknowledge this 
contribution. 

Consider the sequence {sf;} where 5^; = — . It is easy to verify that this sequence is 

positive. Notice that {sk}, where 



« = «+i 



1 

2^ 



2/2^' 



is a positive sequence, but {m *), where ^^ = 1, ^ e Mo, is not a positive sequence in general 
for M e IR . Indeed, if we put u--\ then we arrive at a negative sequence. 

It is easy to see that the difficulty pointed out by this example is that u can take nega- 
tive values. Restricting u to nonnegative values in II (i.e. the Stieltjes moment problem) 
would have made the converse of the above theorem to also hold. Therefore, it should not 
be surprising that for even powers of u stronger results hold. Now we give the follow- 
ing corollaries which give further insights on the characteristics of the sequence {(k) for 
extracting a positive subsequence. 

Corollary 5.5. Let [sk] be a positive sequence. Then the sequence [sk] defined by s^ — s^+t 
for any ( € 2Mo w cf positive sequence. 

Proof. Note that the constant sequence {u ] for ( e 2Mo is positive for any m e M. Apply 
Theorem l5.4l n 



Corollary 5.6. If the sequence [u *) is positive then 

(1) tie 2Fo /or every i e 2Fo, and 

(2) for every odd positive integer i, l\ 



A-i+4 



+1 
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Proof. By Corollary 15. 51 since every tail of the sequence {m *) with {k starting from an even 
integer is positive, 

u^' > for every m e R. 

Hence tj e 21No for every / e 2Fo. 

Again using Corollary 15 .5 1 for every (/ - 1) e 21No, the two by two matrix of the form 






/- 

/.+i 



u ' u 
must have nonnegative determinant. Then 

(26) /'-i+^'+i-M^^' >0. 

Since the inequality (|26T l is true for every m e R, the statement (2) follows. n 

The converse of this corollary does not hold in general. To see this, consider the se- 
quence {m*) W\i\\{(k\ - {0,1,2,4,6,...). The determinant of the corresponding 3 x3 Hankel 
is -u (u - l)", which is a negative function. 

Theorem 5.7. Let {sk} be a positive sequence. The subsequence {sk] given as Sk — Sk+i^ is 
positive if (k — kd + (ofor any d € Mo, and (q e 21No. 

Proof. Suppose i^ - kd + {q where d e INq, and ia e 2]No. Then for all real numbers 

XQ,x\,X2,- ■ ■ ,x,„, we have 



/ J XiXjSi+j - 2_^ XiXjS 



i+j+(i+j)d+eo 



i,j=0 



i,j=0 



m „ 

= Y, XiXj u<'^''^'^<^'^''^J^'odcr(u) 

i,j=Q -'» 

r \ { '" 



n+d)i 



do-{u) > 0. 



The above results show that any subsequence of a positive sequence is positive if they 
are extracted in a certain periodic manner A converse to these results will be proved 
in a later section using matrix completion method. Specifically, we will establish that 
extraction of a subsequence from a positive sequence must be in a periodic manner for the 
subsequence to be positive. Here we provide an example. Consider the sequence given by 

1 

This sequence is positive since, for every m, the determinant of the Hankel matrix given 
by 



A„ = 



1 
e 
1 

2? 



1 



1 
2? 

1 
3? 



1 



1 



(OT+l)e'"+' 
1 



(m + 2)e 



1 



.m+2 



(OT+l)e'"+i (m + 2)e 



m+2 



(2m+\)e 



2m+l 



p(m+l)2 



n 

!=0 



(2/ +1)1 



2/ 
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is positive. But the subsequence given as 

1 

is not positive since the terms are not taken in a periodic manner as defined by the above 

theorem. We can easily verify this by looking at the determinant of the second Hankel 

matrix: 

1 1 

e Ip^ 

02= <0. 

1 1 

6. Sub-Moment Problems 

Given any two positive sequences \s]^] and \~Sj^ with moment solutions cr{u) and &{u) 
respectively, o-(u) - &{u) if and only if Sk - h for all k. Therefore the moment solutions to 
any two distinct positive sequences can never be equal. If {5jt) is a positive subsequence of 
{sk} with corresponding moment solutions d-(u) and cr(u) respectively, then these solutions 
can not be the same. In this section, we develop techniques to compare a-(u) and a-(u). 

Consider a positive sequence {sk} with its moment solution cr(u). If some of the terms 
from {sk} are missing then we obtain a subsequence which, if positive, gives another mo- 
ment solution. For an appropriately chosen {£k] £ INo, assume that {s/c} = {i*:+4) is positive 
and has the moment solution d-(u). Write 

(27) r u^'do-iu) = I u''^^^do-(u). 

Jr Jr 

n 

Then for any polynomial P(u) - y a^u ,we have 

/t=0 

(28) \ P(u)dd-(u)^ f Pe,(u)do-(u), 

JR JK 

n 

where P4(m) = ^fli;M*^+^^ 

k=0 

Here, we present some results concerning moment solution to a shifted sequence Sk - 
Sk+i,€ e 2E^o, which has already been proved to be a positive subsequence of a positive 
sequence Sk- 

Theorem 6.1. Let [sii] be a positive sequence and [sk] be its subsequence given as s^ — 
Sk+t for a fixed { e 2E^o- Let criu) and &{u) respectively be the moment solutions of the 
sequences {sk] and {sk}. If one of the moment problems is determinate, then 

(29) I f(uWdo-(u)^ I f(u)dd-(u), feLlnll. 
Jr Jr 

Proof To prove the first assertion, define the following two linear functionals: 

<1>i(/(m))= f /(m)mW(m), <l)2(g(M)) = I g{u)dcr{u) 
Jr Jr 

for f ^LTg- and g e L^^.. To see Oi and ^2 are bounded, we have 



<Dl (/(«))< f\f(u)uycr(u 
Jr 



) < \\f(u)\\r2 ^i^t<^ 
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and 



3)2(/(«)) < r \g{u)\d&{u) < \\g{u)\y v^ < . 



Now for any function f e L^-H L? , 

<!>(/(«)) = 0i(/(m))-02(/(«)) 

is a bounded linear functional. 

Observe that for any polynomial P(u) = a„,M™ H h oim + aq, 

(D(P(m))= I P(u)u^do-(u)- I P(u)dcr{u) 



u)- I P(m) 
Jb. 

I \a„u'"^ H Va\u^ +aQU \dcr(u)- I [a„u'^ -\ h a i m + ao] '^5"(m) 

Jr Jr 



=am(Sm+l - Sm) + ■ ■ ■ + a\{si+[ - s\) + aoisc - So) 
=0. 



By Corollary 14.11 and due to the determinacy condition of the moment problems, the set 
of polynomials is dense in L^ or L? . Hence the set of polynomials is dense in L"^ n Lj. 
Therefore, by the Hahn-Banach Theorem, O is identically zero. Thus equation (|29] l holds. 

D 

The above theorem is particularly useful when the first finitely many moments are miss- 
ing. It says that we can find the moment solution with the remaining data and reconstruct 
the missing moments by using formula ( |29] l. Whether the moment solution of a subse- 
quence is determinate or not is to be studied more carefully in comparison with the deter- 
minacy of the moment solution of the original sequence. It is, however, known that the first 
term in the subsequence [skf] can be modified to make the moment problem determinate, 
but the solution will not remain the same. 

As a remark to Lr space of sub-moment solutions, the following theorem gives us a rela- 
tion between Lr spaces related to a moment sequence and its sub-moment sequence. Here 
we require the evenness assumption on I to guarantee that v is a nondecreasing measure. 

Corollary 6.2. Consider a positive sequence {sk] and its subsequence {S^.) = {si^+c} for 

{ 6 2]No. Let cr(u) and d'(u) be the moment solutions of the sequence and the subsequence 

respectively. Then there is a measure v{u) absolutely continuous with respect to cr(u) such 

that 

(30) (4n4)cL2. 

Proof. Let f eL^nLl.. By equation ( |29] ) we have 

I f(u)u^do-(u)- I f(u)dd-(u), feLl-HLl.. 
JR JR 



Set dv - u dcr. Then v << cr and 



iu)dv(u)= I f(u) 
Jn 



(31) f(u)dv(u)= I f(u)dcr(u). 

Jn 

Hence /eLy. n 

Now consider a general sub-moment sequence {sk} given as Sk - s<r+4 for an appropriate 
{{k) £ INo- It was shown in the previous section that for any k, 

(k = kd + to, where def^o and /"o, /"i e 21No. 
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Modifying the functional <t)i in the proof of Theorem l6. 1 I gives the following result. 

Theorem 6.3. Let [s^] be a positive sequence and [s^] be its subsequence given by s^ — 
Skd+Co- Let cr(u) and (5"(m) respectively be the moment solutions of the sequences [sk] and 
[sk]- If one of the moment problems is determinate, and f(u ) e L^for any f e L'^, then 

(32) r f(uV''do-(u) = r f(u)dd-(u), feLlnll. 



R 



A precise relation between L spaces of moment solutions to the original sequence and 
that of one of its positive subsequences can be useful in characterizing the sub-moment 
solutions and hence approximating and/or reconstructing the missing data. 

Exploring equation ( [32] i a little further, for A € C,^ =3/1^0, define 

f(u) = -^. 



Then we have 



\f{u)\^dcr(u)^ ■ — < — = — <oo. 



since \u- A\> \y\. Therefore, / e L^. By the same argument, / e L? . Similarly, it can be 
shown that f(u'') e L^. Then by TheoremlOl 

, u'^°do-(u) f dd-(u) 



r u'^oda-(u) _ r I 



u — A 

Applying Theorem 14.41 to the equation (|33T l yields the following result which connects a 
moment problem with the polynomials corresponding to its sub-moment problem. 

Theorem 6.4. Let [sk] be a positive sequence and [sk] be its subsequence given as sk — 
Skd+iQ- Let <t(u) and d'iu) respectively be the moment solutions of the sequences [sk] and 
[sk]. If the moment problem of{sk] is indeterminate, then 



(34) 



r u'^o 

Jr w 



u^^dcriu) _ A{A)(I)(A) - C(A) 
B{A)(p{A) - D(Ay 



where A(A),B(A),C(A),D(A) fonn a Nevanlinna matrix of the sub-moment problem, and 
(heN. 



Next, for the shifted sub-sequence we discussed in Theorem l6.1l we want to investigate 
determinacy of its moment problem of a special sub-moment sequence in relation to deter- 
minacy of the original moment problem. Recall that the limit circles Koa{A),A € C, 3/1 ^ 
defined in the Section 2 provide a fundamental concept for studying determinacy of mo- 
ment problems. 



Theorem 6.5. For the same moment sequences and conditions as in Theorem 16.71 for 
/I € C, 3/1 ?i the following holds. 

(35) rdoM^^^^, rdoM 

JR u-A Jr. u-^ 

where C is a constant depending on A. 

Proof. It is already shown that, for 

/(m) = ^— , /leC,3/l?iO, 
u — A 
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/ e 4 n 4. Then by TheoremO 



u dcr(u) r dd-(u) 
u — A Jr u- a 

But = R(u) + , where R(u) = V A'u^~^~'. Therefore, 

u- A u — A ^ 

1=0 

C d&Ou) C , , , , , f r dcr(u) 
— —^ R{u)d&(u) + r — —. 

JR U-Z Jr Jr U-A 



Hence 



R 



d&iu) 
u — A 



= C + A^ 



R 



dcr(u) 
u — A ' 



where C = ^^' 
Let 



Si-l-i- 



i=0 



i) 



f dcr(u, 

Wct{A) = and wg-iA) 

JR u-^ 

Then Theorem l6.5l states that 

(36) W5-(/l) = C + /w^(i), 



r dd-(u) 
jR u-'^ 



for Sk and Sk specified in that theorem. 

Recall that points Wa- and Wg- lie on the circumferences of the circles Koo{A) andA'oo(/l) 
corresponding to the moment problems of [sk] and [sk] respectively. Equation ( |36] | pre- 
cisely describes the distortion of the circle corresponding to the original moment sequence 
in relation to its tail. 

7. Completion of Positive Hankel Matrix 



Recall from Definition l2.1l that the positivity of an infinite sequence {sk) is equivalent to 
the positivity of Hankel matrices of the form 



H„ 



so 
s\ 



■Sm+1 



Sfii 
Sm+ 1 

S2m 



Determining whether a subsequence of a positive sequence {s^} is positive is then equiv- 
alent to a matrix completion problem. A partial Hankel matrix is a partial matrix that is 
Hankel to the extent to which it is specified. All the specified entries lie along certain skew- 
diagonals (positive sloping diagonals), and all the entries along a specified skew-diagonal 
have the same value. For example, for a moment sequence {s^} with every (2; -1- l)-th entry 
missing we have a corresponding Hankel matrix 



^5 



^0 


* 


S2 


* 


S4 


* 


* 


^2 


* 


54 


* 


^6 


^2 


* 


S4 


* 


S6 


* 


* 


54 


* 


^6 


* 


^8 


S4 


* 


^6 


* 


Si 


* 



S6 



S& 



SlO 
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Essentially, we want to determine whether a particular partial Hankel matrix can be com- 
pleted to a positive Hankel matrix. Note that here we are interested in completing the 
partial Hankel matrices corresponding to a sparse moment sequence. 

The problem of matrix completion has been studied in great details (see ^, O, |l6l, 
II20I . for example). Most common matrices that are of interest in matrix completion are 
positive definite matrices, positive semi-definite matrices, co-positive matrices, totally pos- 
itive matrices and Toeplitz matrices. Matrix completion properties have several applica- 
tions in sparse moment problem. Some results on multidimensional trigonometric moment 
problem have been proved using moment completion in [2 1. 

Most of the results in matrix completion deal with the case when the main diagonal 
entries are specified. This creates a problem in completing a partial Hankel matrix. When- 
ever a (2/)-th term of a moment sequence is missing, the corresponding Hankel matrix has 
a diagonal entry missing. 

Definition 7.1. Given an nxn matrix A, let w,v be subsets o/lNo. A submatrix A[w,v] 
consists of the entries in rows w and columns v. The principle submatrix is A[w] — A[w,w'\. 

For example, 

12 3 4 



5 6 7 8 
9 10 11 12 
13 14 15 16 



A[{2,4)]=A[{2,4),{2,4)] = 



6 8 
14 16 



Matrix completion encounters an important inheritance structure. If A is a matrix with 
property X, then in most cases of matrix completion of interest, the fully specified principle 
submatrices of A inherit the property X. Positive (semi) definite matrix completions have 
this inheritance structure. But the converse is not true for a general positive (semi) definite 
matrix. That is, a partial matrix with all of its fully specified principle submatrices positive 
(semi) definite does not necessarily have a positive (semi) definite completion, |20|. Simi- 
larly, if a partial matrix A admits a completion of rank < k, then every specified submatrix 
of A has rank < k. 

With the inheritance structure of matrix completion, the matrix completion problem for 
positive Hankel matrices becomes equivalent to finding all the positive subsequences of a 
moment sequence. That is, the principle sub-Hankel matrices of a positive Hankel matrix 
are those corresponding to positive subsequences of a moment sequence. 

Given a Hankel matrix H corresponding to a moment sequence [sk] consider a fully 
specified principle submatrix H[a\, where a - {aQ,a\,a2,...,an] £ INq. For H{a] to be 
also a Hankel matrix we require that 

a,-, + aji = a/j + "jj whenever ii + ii = ji + jj. 

Then for any k, we have 

ak-i + ak+i 



2 

Also note that every principal submatrix H{a\ has (1,1) entry S2i for some /. With this 
observation and Theorem 15. 71 we have proved the following result: 

Tlieorem 7.2. Let [sk] be a positive sequence. A subsequence [sk] given as Sk — Sk+ii^, is 
positive if and only if£k — kd + £ofor some d e INo ond {q e 2]No. 

In terms of matrix completion the above theorem gives a necessary and sufficient condi- 
tion for completion of a partial Hankel matrix. Interestingly, this result for Hankel matrices 
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seems to be similar with the result for Toeplitz matrices proved in fJTl. The reader is re- 
ferred to |3 1 for more on matrix completion problem for Hankel and Toeplitz matrices. 

8. Remarks 

The problem of identifying positive subsequences of moment sequences and the cor- 
responding matrix completion problem is the forward problem associated with sparse se- 
quences and sparse matrices. The corresponding inverse problem of extension of positive 
sequences as sparse subsets of moment sequences is closely related. The necessary and 
sufficient conditions obtained suggests that positive sparse sequences can be completed 
by extending the measures associated with these sequences by Theorem |6.4| and calculat- 
ing the arising moments. This suggests that the moment sequence can be recovered to a 
great degree from its positive sparse subsequences. The entropy concept introduced by J. 
Chover f8l for extensions of positive definite functions was used to find a unique extension 
of positive definite functions which maximizes the entropy. In parallel to this classical and 
well-known theory, we anticipate that for each positive sequence arising from the Nevan- 
linna class there is a unique extension in the Nevanlinna class that maximizes the entropy. 
The corresponding phrasing of this problem for positive Hankel matrices gives a norm, in- 
duced by the entropy, which allows recovery of these structured matrices from their sparse 
positive submatrices. Relationship of this approach to finding the matrix with minimum 
nuclear norm that fits the data in Candes and Recht's |i6J is not clear to us at the present 
time, but is a subject of further investigation. 
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